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' We reexamine the Yang- Yang- Takahashi method of deriving the thermodynamic Bethe ansatz 

, equations which describe strongly correlated electron systems of fundamental physical interest, such 

' as the Hubbard, s — d exchange (Kondo) and Anderson models. It is shown that these equations 

contain some additional terms which may play an important role in the physics of the systems. 
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It is well known that many one-dimensional (ID) and effective ID models describing strongly correlated electron 
systems of fundamental physical interest are diagonalized exactly by following set of the Bethe ansatz (BA) equations: 
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a=l 



exp {ikjL) = Y[ si{uj — Aq) (la) 



ei(Aa - Uj) = - ]^ e2(Acj - A^). (lb) 
i=i /3=1 



Here, e„(x) — {x + in/2)/{x — in/2), kj are electron momenta, Uj = u{kj), N is the total number of electrons on an 
r , interval of size L and M is the number of electrons with spin "down" . The eigenenergy E and z component of the 
total spin of the system are given by 
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iJ^^c.,, S^^-N-M. (Ic) 



^ ' where ujj — oj{kj) are electron energies. 

00 At LL!{k) ^ k'^, u{k) ^ k and uj{k) ^ — cosk, u{k) ^ sinfc, Eqs. (1) correspond to a ID electron gas with (5-function 

interaction ||l|,|| and the ID Hubbard model respectively. In the theory of dilute magnetic alloys |^-p, these 
^ ] equations describe the excitation spectrum of a free host in terms of interacting Bethe excitations, while impurity 
■ terms, omitted in Eqs. (1), account for an elastic scattering of Bethe excitations at the impurity site. The case of 
^\ Lo{k) = fc, u{k) = corresponds to the s — d exchange (Kondo) model while the case of uj{k) = k, u{k) fc^ 

0^ ■ corresponds to the Anderson model j^. In addition, Eqs. (1) can be used to study the Anderson and Kondo impurities 
embedded in a host with a nonmetallic behavior of the density of band states |l^ . 

The equations (1) correspond to a finite system of electrons. A general method of deriving the thermodynamic 
BA equations has been developed by Yang and Yang |l^ for the case of spinless particles, and then generalized by 
Takahashi ||l^ to the case of particles with internal degrees of freedom. 
(~| In this Letter, we reexamine the Yang- Yang- Takahashi method and show that the thermodynamic Bethe ansatz 

O • equations contain some additional terms which may play an important role in the physics of strongly correlated 
electron systems. These additional terms result from a more accurate computation of a contribution of spin degrees 
of freedom to a variation of the system's energy. 

The structure of the BA equations shows that it is natural to treat the system of — M electrons with spin "up" 
and M electrons with spin "down" in terms of N charge excitations with "rapidities" Uj = u{kj), j = 1, . . . , iV and 
' spin excitations (spin waves) with rapidities Aq, a = 1, . . . , Af. In accordance with the Yang- Yang- Takahashi method, 
to derive the thermodynamic BA equations, one needs to note first that in the limit of a sufficiently large physical 
system Eqs. (1) admit bound spin complexes in which spin rapidities Aq are grouped into "strings" of order n in the 
complex A plane, 

A^^) =A2 + ^(n+l-2j), J = !,...,«. 

In addition, for many of the above-mentioned physical systems the BA equations admit also charge complexes in 
which charge excitations are bond to spin complexes |^,^,|l2|. Nevertheless, to avoid more tedious expressions we start 
our analysis with the case in which charge complexes are absent. 
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Then, Eqs. (1) take the form of equations for charge excitation (particle) momenta kj and rapidities of spin 
complexes, A^, a = 1, . . . , M„, where M„ is the number of complexes of order n and hence J2n ^-^n = M, 

oo M„ 

2nNj = kjL + J2 T^PnHkj) - AS], (2a) 

n=l a=l 

N oo Mm 

2nJ: = J2 W - -EE ®r^m{Xl ' Xf), (2b) 

j=l m=l /3=1 

where p„ (x) = tt + 6'„ (x) , and iVj and J" are sets of quantum numbers of the system corresponding to particles and 
spin complexes, respectively. Here, 9n{x) = 2 arctan (2a;/n) and 

min(n,m) — 1 

Qnm.{x) = (1 - S„m)0\ 

n—m\ 

{x) + 9n+m{x) + 2 ^ ^1 

n— m|+2fc 

fe=l 

In the continuous limit, L oo, A'' ^ oo, M ^ oo, but N/L and M/L are constant, Eqs. (2) take the form 
of integral equations for the "particle" [p{k) and cr„(A)] and "hole" [p{k) and o-„(A)] density distributions of charge 
excitations and spin complexes, respectively, 

p{k)+p{k) = — +u'{k)y2 dXan[u{k)-X\an{X), (3a) 



a„(A)= f dkan[X-u{kMk)-f2 f dX' A^miX - X')a^{X'), 

m=l 



(3b) 



where a7i(a;) = (2n/7r)(n + 4a; ) , = du/dk, and the matrix Anm{x) is given by 

min(n,m) — 1 

(x) = SnmS{x) + (1 {x) + 2 a|„_„|+2fe(a;). 

In terms of the densities, the energy of the system is found to be 

= j dkuj{k)p{k) (4) 

The number of unknown functions in Eqs. (3) is twice bigger the number of equations. Needed additional equations 
for fundamental (renormalized) energies of Bethe excitations 

pyk) cr„(A) 

are derived from the condition of a minimum of the thermodynamic potential of the system 

n = E-TS- HS^ - AN, (5) 

where T, S, A, and H are the temperature, entropy, chemical potential, and an external magnetic field, respectively. 
The condition S^l = 0, where the energy's variation is given by 

^^E = j dkw{k)Sp{k), (6) 



results in 



1 f 

e{k) = uj{k) -T^H-A-Y^ dXanHk) - A] F[-Kn{X)\ (7a) 
F[nn{X)]^nH + j dku {k) an[X - u{k)]F[~e{k)] 

OO „ 

+ Y j dX' A„„(A - A') F[-Km{X')] (7b) 
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where F[f{x)] = T\n{l + exp [f{x)/T]}. 

The analysis of the thermodynamic BA equations (7) in the Umiting case of the Kondo system shows that the 
expression (6) is incomplete. Although, in accordance with Eq. (3a) the variation of the density distribution of 
particles is related to the variation of the density distributions of spin complexes by the relation 



Sp{k) + Sp{k) = u'{k) / dAa„[u(fc) - A](5cr„(A), 



the latter contains only the derivative of the charge rapidity u'{k). It is easy to see that if one sets u'{k) = 0, as it 
takes place in the Kondo model, the dependence of the energy variation on the spin degrees of freedom completely 
disappears from Eq. (6), and Eq. (7b) immediately leads, therefore, to unphysical result: All k„(A) are positive, and 
hence spin excitations are absent in the ground state of the system. The expression (6) does not account thus for a 
total contribution of spin degrees of freedom to a variation of the energy of the system. 
In the continuous limit for spin excitations only, Eq. (2a) for a particle momentum, 

27r oo „ 

—N, = kj + J2 d\pn[u{kj) - A]a„(A), (8) 

n=l •' 

clearly shows, however, that the momentum of a charge excitation, and hence its energy to^kj), are determined by 
both the quantum number of charge excitations Nj and the density distributions of spin complexes ct„(A) even if 
u{k) = 0. Therefore, the total variation of the system energy should be written in the form 

y6E^ I dkuj{k)6p{k) + I dkp{k)SLj{k), (9a) 



L 

where SLL}{k) is the variation of the energy of a particle at fixed quantum numbers of charge excitations, SNj = 0, 

du;(fc) dujjk) ^„/dAp„[u(fc) - A](5g„(A) 

'^^ ' dk dk l + 27ru'(fc)^„/dAa„[A-u(fc)](T„(A)' ^ ' 

Making use of Eq. (3a) for the denominator in the right-hand side of Eq. (9b) , one easily finds the final expression 
for the total variation of the energy of the system. 

Thus, while the energy of the system is given by Eq. (4a), its variation contains two terms. The first term corresponds 
to a contribution of a variation of charge quantum numbers at fixed spin quantum numbers of the system, while the 
second one corresponds to a contribution of a variation of spin quantum numbers at fixed charge quantum numbers. 
Taking into account Eq. (10), we immediately find that the term 

dk duj{k)/dk 
27r 1 + exp [e[k)/ 1 \ 

should be added to the right-hand side of Eq. (7a). Then, in the Kondo limit, where u{k) — 0, du){k)/dk ~ 1, and 

dk 1 N 
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27r 1 -Hexp[e(A:)/r] L 



Eq. (7b) with the extra term (11) correctly reduces to the standard equation describing the thermodynamics of the 
Kondo model |J|. 

If the Bethe spectrum of a systems contains also charge complexes, Eq. (8) takes the form 

2 oo „ 

—N, ^k,+Y dXp4u{k,) ~ A] [a„(A) + a'^iX)] , (12a) 

where cr^i(A) stands for the density distribution of charge complexes of order n. The total variation of the energy of 
the system is then found to be 
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I dku:{k)5p{k)- I g—^^l^j^f; I d\p4u{k)-X]S[a^iX)+aUX)]- (12b) 

The Yang- Yang- Takahashi method results then in the foUowing set of the thermodynamic Bethe ansatz equations: 

1 °° f 

e{k) = co(fc) - H-A-J2 dXanHk) - A] {F[-k„(A)] - FK„(A)]} , (13a) 
F[n4X)]^nH-J g--^^^|^p„Kfc) - A] + J dku' (k) a^X - u{k)] F[-e{k)] 



OO 



m— 1 



V dX'A,,„,{X-X')F[-K^{X')], (13b) 



^^Kn(A)]=enA)-2nA- / ^ J^^WM 

7 27r 1 + exp [^(/cj/TJ 

/OO „ 
d/cu'(ft)a„[A - M(A:)]F[-e(fc)] + J2 dXA„^{X - A')i^[-U(A')]- (13c) 
m— 1 

where Cn(A) = T^ln [''■n(A)/cr,'j(A)] are the renormahzed energies of charge complexes, while <^i*'''(A) are their bare 
energies, which are different for different physical systems. 

In summary, the additional terms in the thermodynamic BA equations, derived in this Letter, should essentially 
affect the thermodynamic properties of a system at finite temperature. It is clear that the additional terms affect 
essentially also the low-energy physics of such a system, in which particles (charge excitations unbuilt in charge 
complexes) essentially contribute to the ground state of the system, as it takes place in the Kondo and Hubbard 
models. If unpaired charge excitations are absent in the ground state of a system, as it takes place in the Anderson 
model ^,^,|l3| in the absence of an external magnetic field, the additional terms disappear from the zero-temperature 
limit of Eqs. (13). In a weak magnetic field, the ground state of the Anderson model contains also a small portion 
of unpaired charge excitations, and the additional term in Eq. (13c) should result in some small corrections to the 
standard solution. 
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